New Universality Class in the S=l/2 Fibonacci Heisenberg Chains 
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Low energy properties of the 5 = 1/2 antiferromagnetic Heisenberg chains with Fibonacci ex- 
change modulation are studied using the real space renormalization group method for strong ex- 
change modulation. Using the analytical solution of the recursion equation, the true asymptotic 
behavoir is revealed, which was veiled by the finite size effect in the previous numerical works. It is 
found that the ground state of this model belongs to a new universality class with logarithmically 
divergent dynamical exponent which is neither like Fibonacci XY chains nor like XY chains with 
relevant aperiodicity. 
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The magnetism of quasiperiodic systems has been the 
subject of continual studies since the discovery of qua- 
sicrystals in 1984[1]. This problem has been attract- 
ing renewed interest after the synthesis of magnetic qua- 
sicrystals with well- localized magnetic moments [2]. The 
artificial formation of one and two dimensional quasiperi- 
odic structure is also coming into the scope of exper- 
imental physics[3, 4] thanks to the recent progress of 
nanotechnology and surface engineering. Possiblly mo- 
tivated by these experimental progress, the theoretical 
investigation of the quantum magnetism in one and two 
dimensional quasiperiodic systems are started by many 
authors [5-10]. 

The 5=1/2 Fibonacci XY chain, which is mapped 
onto the free fermion chain, has been studied extensively 
by Kohmoto and coworkers [11] by means of the exact 
renormalization group method from the early days of 
quasicrystal physics. It is shown that the ground state of 
the XY chain with Fibonacci exchange modulation is crit- 
ical with finite non-universal dynamical exponents. This 
approach was extended to include other types of aperiod- 
icity and anisotropy[12]. It is clarified that the criticality 
of the Fibonacci XY chain emerges from the marginal 
nature of the Fibonacci and other precious mean aperi- 
odicity in this model. For the relevant aperiodicity, more 
singular behavior with a divergent dynamical exponent 
is realized even for the XY chain [12]. 

On the other hand, the investigation of the S = 1/2 Fi- 
bonacci Heisenberg chains started only in late 90's. The 
ground state of the uniform S — 1/2 Heisenberg chain 
is exactly solved by the famous Bethe ansatz method[13] 
and is known to be in the Luttinger liquid state with 
conformal invariance. This implies that the dynamical 
exponent z is unity and the specific heat C and suscepti- 
bility x behave as C ~ T and x ~ const, at low temper- 
atures. This exact solution is related to the transfer ma- 
trix of the 2-dimensional classical 8 vertex model which 
can be solved exactly. [13] On the quasiperiodic lattice 
also, some 2-dimensional classical models are known to 



have exact solution. [14] However, the exact solution of 
the Fibonacci Heisenberg chain is not derived from these 
exactly solvable vertex models. Therefore we must re- 
sort to the renormalization group approach to clarify the 
reliable low energy asymptotic behavior. For weak Fi- 
bonacci modulation, Vidal and coworkers[7] have shown 
that the Fibonacci modulation is relevant on the basis of 
the weak modulation renormalization group calculation. 
The present author carried out the DMRG calculation 
and investigated the scaling properties of the low energy 
spectrum[5, 6]. In the present work, we employ the real 
space renormalization group (RSRG) method[15], which 
is valid for the strong modulation, to elucidate the ground 
state properties of the S = 1/2 antiferromagnetic Fi- 
bonacci Heisenberg chains. Surprisingly, the finite size 
scaling formula, which fitted the DMRG data in ref. 
[5] well, turned out to be the artifact of the finite size 
crossover effect. The true asymptotic behavior is first 
revealed by the exact solution of the recursion equation 
obtained in the present paper. It is also explained why 
the DMRG data are well fitted by the formula assumed 
in ref [5] within the appropriate range of the system size. 
Our Hamiltonian is given by, 

N-l 

H=J2 J^SiSi+i, (J ai > 0, a, = A or B), (1) 
t=i 

where Si's are the spin 1/2 operators. The exchange cou- 
plings J ai (= Ja or Jg) follow the Fibonacci sequence 
generated by the substitution rule, 



A -» AB, B -> A. 



(2) 



If one of the couplings Ja or Jg is much larger than the 
other, we can decimate the spins coupled via the stronger 
exchange coupling and calculate the effective interaction 
between the remaining spins by the perturbation method 
with respect to the weaker coupling [15]. This type of dec- 
imation scheme has been used to investigate the magne- 
tization process of the Fibonacci Heisenberg chains [8]. 
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FIG. 1: The decimation procedure for J A << Jg- 
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FIG. 2: The RSRG scheme of the Fibonacci Heisenberg chain. 
The letters A and B correspond to the bonds and the up 
and down arrows to the spins which survive decimation. For 
J j± » Jg the leftmost spin and bond in the parenthesis do 
not appear. 



Here we apply this scheme to find the fixed point which 
governs the ground state in the absence of magnetization. 

In the present work, we concentrate on the case of 
strong modulation, max( J A / Jg, Jg/J A ) >> 1. For 
J A << Jg, the spins connected by the Jg-bonds are 
decimated as shown in Fig. 1. The spin- 1/2 degrees of 
freedom survive on the sites in the middle of the sequence 
AA. The two kinds of sequences of bonds are allowed 
between two alive spins, namely ABA and ABABA. Be- 
tween these two alive spins, there exists one singlet pair 
in the former case while two singlet pairs exist in the lat- 
ter case. Therefore the effective coupling is weaker for 
the latter case. This decimation process replaces the se- 
quence ABABA by A' and ABA sandwiched by two As by 
B' resulting in the sequence B'A'B'A'A'B'A'B'A'... Ex- 
cept for B' at the leftmost position, this sequence again 
gives the Fibonacci sequence as schematically shown in 
Fig. 2. As seen from the change of the number of the 
bonds by one step of decimation, this procedure essen- 
tially corresponds to a 3-step deflation. The rigorous 
proof will be reported in a separate paper[16]. 

In the case J A » Jg, the decimation precesses are 
shown in Fig. 3. The three spins connected by succes- 
sive A bonds form a doublet which can be described as 
a single spin with magnitude 1/2. The spins connected 
by the isolated J A bonds are decimated. Therefore this 
decimation process again corresponds to the replacement 
ABABA — » A' and ABA^ B'. In this case, the result- 
ing sequence is the exact Fibonacci sequence. After the 
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FIG. 3: The decimation procedure for J A >> Jg- 



first decimation the A' bond becomes weaker than the 
B' bonds. Therefore the decimation rule for the case 
J A « Jg applies for the further decimation procedure. 

The effective coupling can be calculated by the 
straightforward perturbation theory in weaker coupling. 
For the n-th iteration, we have, 
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for Jg >> J A , (4) 



for J A » Jg, (5) 



where the variables with refer to the values after n- 
step iteration. 

The ratio j!™VJtT^ decreases under renormalization 



as, 
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(n > 1). 



(6) 



This implies that the perturbation approximation be- 
comes even more accurate as the renormalization pro- 
ceeds. Therefore the aperiodicity is relevant in consis- 
tency with the result of the weak modulation renormal- 
ization group method [7]. Taking the both results into 
account, we may safely expect that the ground state of 
the Fibonacci Heisenberg chain is governed by the strong 
modulation fixed point obtained in the present approach 
in the entire parameter range of J A / Jg ^ 1 . 

The solution of the recursion equation (3) is given by, 



J A = J A 
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> for Jg >> J A , (7) 
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The length of the 3n-th Fibonacci sequence is equal to 
the Fibonacci number F^ n which grows as <f> 3n for large 
n where (j) is the golden mean (= 1+ 2 V ^ ). Therefore the 
chain of length N ~ <f> 3n reduces to a single pair of spins 
after n decimation steps. This implies that the smallest 
energy scale AE for the finite Fibonacci chain with length 
N scales as, 



AE < 



exp(-(lnA^/31n< ? !)) 2 ln2) 



-K(lnAf) 2 _ jy- K \llN 



with k = ln2/(31n<£) 2 (9) 
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for large enough TV, irrespective of the value of — — . It 

J B 

should be noted that the dynamical exponent diverges 
logarithmically. 

This size dependence implies that the number of the 
magnetic excited states with energies in the interval 
AE ~ AE+dAE the magnetic excitation ND(AE)dAE 
scales as, 



This low temperature behavior should be contrasted with 
that of the uniform 5=1/2 antiferromagnetic Hcisen- 
berg chain C ~ T and \ ~ const, which is less singu- 
lar than the present Fibonacci case. This is due to the 
logarithmic divergence of the dynamical exponent in the 
present case. To check the reliability of the present RSRG 
scheme, we also applied the same procedure for the Fi- 
bonacci XY chain to find 



ND(AE)dAE 



AE ~ TV" 



(15) 



/ ^cxp (^>aV) J d (^exp [-^m^j 



" X ' I V( ' X1> I ^A^E 



1 



' \ b n Ah) dAE 



(10) 



with a scaling function f(x). Because the density of state 
per site D(AE) should be finite in the thermodynamic 
limit TV — > oo, the scaling function f{x) tends to a finite 
value as x — > oo. Therefore we find, 

D{AE)dAE 



2nAE 



'V: j \lhn^\dAE (11) 



for large enough TV. Accordingly, the low temperature 
magnetic specific heat C should behave as, 



C 



d f T 

TV / AED(AE)dAE ~ NTD(T) 
Jo 
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(12) 



The magnetic susceptibility at temperature T should be 
the Curie contribution from the spins alive at the energy 
scale T. The number n s (T) of such spins is given by 



n s (T) ~ 2 TV / D(AE)dAE, 



(13) 



because two spins are excited by breaking a single effec- 
tive bond with effective exchange energy less than ksT. 
Therefore the low temperature magnetic susceptibility \ 
behaves as, 
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(14) 



2 / (Jk J r \\ 
with z = — — -In max -P~. — — . This reproduces 

the exact result by Kohmoto and coworkers [11] in the 
limit max(-^,-^| >> 1. Therefore our RSRG 

scheme is reliable at least for strong modulation. 

The present results appears to be in contradiction with 
the results of ref. [5] , in which the present author carried 
out the DMRG calculation for the Fibonacci antiferro- 
magnetic Heisenberg chains. In ref. [5], we performed the 
finite size scaling analysis of the lowest energy gap AE 
based on the assumption that it will behave in the same 
way as the XY chain with relevant aperiodicity, namely 
as AE ~ exp(— cTV"). However, the present analysis sug- 
gests the different behavior. Although we tried to replot 
the previous data using the scaling (9), the fit turned out 
to be very poor. The reason of this discrepancy will be 
understood in the following way. 

For finite J^/Jg, the perturbation approximation re- 
quires the higher order corrections which modify the re- 
cursion equation (3) in the form, 

j(n)3 

t(t»+1) _ A f , ( j{n) I r(n)\ 

J A - "7H2 A /J B >' 



J B 



'B 

T (n)2 



—MJ A /J B )■ 



(16) 
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It should be noted that the correction factors /a and /b 
depend only on the ratio J^/Jj^ and satisfy /a(0) = 
1/4, /b(0) = 1/2. This leads to the recursion equation 
for XW - (X<g\ X§>) = (lnJ^lnjj^) as 

= ( 2 -1 ) X(n) + rt X A ~ X B } ) 

with (j, — (/ia,Mb) = (hi/A, hi/e)- If the function n(X) 
is approximated by a linear function of X as fi(X) — 
jX + no with 7 = (7a, 7b), we have 

X<" +1 ) = M m X^+no, (18) 
where M m is a 2 x 2 matrix 

u„ = ( i m -?t 7A V ( i9 ) 

2 - 7s -1 +7s / 
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One of the eigenvalue of M m is unity. If another eigen- 
value A m (= 1 + 7b — 7a) is larger than unity, the solution 
of (18) grows with n as 

X {n) oc \ n m . (20) 

In this case, both hijj^ and hiJg^ scale as A^. There- 
fore the lowest energy scale of the chain of length TV also 
scales as, 

AE ~ exp(-C*A|T*) ~exp(-CWW) 

~ exp(-CN") with u = (21) 

within appropriate range of system size N. This is the 
reason why the behavior (21) is observed in DMRG cal- 
culation for finite systems. We have numerically diag- 
onalized the Hamiltonian of the clusters BAABABAAB 
and BABAABAABAB which reduce to a single A'-bond 
and B'-bond after decimating B bonds. Using these nu- 
merical data, it is verified that the effective value of 
X m is larger than unity although it actually depend on 
J^/Jg. As the renormalization proceeds, of course, the 
ratio Jp^/J-Q decreases and the true asymptotic behav- 
ior (9) is reached. More details of this calculation arc 
presented in [16]. 

This crossover behavior implies that the extremely low 
temperature is required to observe the true asymptotic 
behavior (12) and (14) in weak modulation regime. In- 
stead, the behaviors expected from (21), namely C <~ 
l/(lnT) 1+1 / w and X ~ l/(T(lnT) 1 / w )[12], would be ob- 
served in the intermediate temperature regime. 

In summary, using the RSRG method, we have shown 
that the ground state of the 5 = 1/2 Fibonacci Heisen- 
berg chain belongs to a new universality class in which 
the energy gap scales as exp(— «(ln7V) 2 ) where k is a uni- 
versal constant independent of modulation strength. The 
low temperature behavior of the magnetic specific heat 
and magnetic susceptibility is predicted. The relation- 
ship to the previous numerical results [5] which appear to 
contradict with the present calculation is also discussed. 
The details of the calculation and proof will be reported 
in a separate paper, which will also includes the discus- 
sion of the general XXZ case[6] and ground state phase 
diagram [16]. 

We have found a new quantum dynamical critical be- 
havior (9) which was so far unknown in the field of quan- 
tum many body problem. Similar 'singular dynamic 
scaling' is, however, known since 80's for the classical 
Ising model on the percolation clusters with Glauber 
dynamics. [17] In spite of the geometrical self-similarity 
common to this classical model and our quantum model, 
they look very different in many aspects. Although the 
underlying physics is still unclear, further investigation 
on this point might lead to a more profound understand- 
ing of both systems. 



After this work is completed, the preprint by Vieira[18] 
appeared in e-print archive in which some of the present 
results are derived. In addition, Vieira satisfactorily ap- 
plied this method to the Heiscnbcrg chains with relevant 
aperiodicity. Similar approach is also applied to the two 
dimensional quasicrystal.[10] We thus expect the RSRG 
method is widely applicable to various problems in the 
field of quantum magnetism in quasiperiodic systems. 
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